Criticality and oscillatory behavior in non-Markovian Contact Process 
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A Non-Markovian generalization of one- dimensional Contact Process (CP) is being introduced in 
which every particle has an age and will be annihilated at its maximum age r. There is an absorbing 
state phase transition which is controlled by this parameter. The model can demonstrate oscillatory 
behavior in its approach to the stationary state. These oscillations are also present in the mean- 
field approximation which is a first-order differential equation with time-delay. Studying dynamical 
critical exponents suggests that the model belongs to the DP universlity class. 
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I. INTRODUCTION 

Studying phase transitions in the systems far from 
equilibrium has been a topic of growing interest in re- 
cent years [Qj2| . Specially systems with absorbing states 
which can not evolve further once they are trapped in one 
such state have been an interesting subject of research. 
Various models with one or more absorbing states have 
been studied which belong to a few universality classes 
and mostly to that of Directed Percolation (DP). Accord- 
ing to the DP conjecture |j|,f§], every phase transition in 
a system with a single absorbing state having short range 
interactions with no special symmetry or quenched dis- 
order H belongs to the DP class. There has been exam- 
ples of absorbing state phase transitions without some 
of the DP conjecture conditions that still belong to the 
DP class such as systems with an infinite number of ab- 
sorbing states IK?]] . There has also been another uni- 
versality classes for parity conserving systems and as 
recently proposed, for systems with infinitely many ab- 
sorbing states coupled to a conserved field Jj|[l0| . 

Although the Markovian property have been implicitly 
accepted in these models it is not essential for a nonequi- 
librium phase transition. Usually, adding some kind of 
memory to the system, such that the system should refer 
to its history in order to define its future, gives rise to 
some new interesting behaviors that are absent in Marko- 
vian ones |Il],[l2|]. However properties of non-Markovian 
nonequilibrium phase transition have not been studied. 

In this paper a non-Markovian variant of the Contact 
Process (CP) O] is introduced and the critical behavior 
is investigated. Standard CP in its continuous time ver- 
sion is a lattice model in which every empty site is being 
occupied by a particle with rate \n/z and every parti- 
cle is removed with rate one, where z is the coordination 
number, n number of occupied nearest neighbors and A 
a positive parameter controlling the creation rate. The 
system has a second-order critical point at A c = 3.2978 
and it belongs to the DP class 

In this model I introduce a memory for each particle: 
every particle knows when it has been created. Like stan- 
dard CP, every site is being occupied by a particle at a 



rate proportional to the number of its occupied nearest 
neighbors, while every existing particle will die exactly 
at age r. For large values of r the particles live enough 
to reproduce plenty of another ones and the system can 
remain in its active state. As r is decreased, each particle 
has less time to create another particles and for r < t c 
the system will be trapped in its absorbing state with 
probability 1, where there is no existing particle and no 
new particles can be born. 

Attributing age to the particles in the CP model has 
been suggested earlier but not in a way that leads 
to a non-Markovian model. Although some interesting 
alterations in the dynamical behavior of the system have 
been observed. 

Non-Markovian property gives rise to some oscillations 
in the density of particles. These oscillations are also sup- 
ported by the Mean-Filed approximation. Because of the 
non-Markovian property of the system, there is a delay 
parameter in the mean field equation and this reproduces 
the oscillatory solutions observed in the simulations. By 
mean-field approach, existence of the phase transition is 
justified and a critical age can be found, but like standard 
CP the critical behavior is not described completely. 

In this paper, the critical behavior of the model is stud- 
ied using the time-dependent Monte Carlo method. The 
critical dynamical exponents has been calculated, and 
shown to be in good agreement of those of DP. 



II. MODEL 

The model is defined on a 1-dimensional lattice and 
with continuous time. Every site is either empty or oc- 
cupied by a single particle. There is a chance for a vacant 
site to be occupied provided there are occupied sites in 
its nearest neighborhood. A new particle is born in an 
empty site with rate n/2 (n number of occupied nearest 
neighbors). Every particle will die exactly at time r after 
its birth. 

Obviously there must be a phase transition in the sys- 
tem with density of the particles as the order parameter. 
For small values of r particles die fast and eventually the 



1 



system is trapped in its absorbing state. For large r's, 
particles have a large lifetime and reproduce sufficiently 
other ones to keep the system active. Fig. 1 shows a 
single cluster in a realization of the model for t = 3.5 
and up to t = 100. As can be seen, every particle has a 
definite lifetime and here the system is in its active state. 



FIG. 1. A typical space-time cluster started with a single 
particle in the origin for r = 3.5 and up to t = 100. 

Like CP, here the creation process is not history- 
dependent, however the death process is, and thus we 
have a non-Markovian process. To find out that how 
many particles are removed at each moment, we should 
know how many particles had been created at time r ear- 
lier. Therefore a knowledge of only the present state of 
the system is not enough to find out the future states. 



III. MEAN-FIELD EQUATION 

Let (Ti (x) and a (x) denote the state of site x, and p be 
the density of particles. is 1 if the site is occupied 

and if it is vacant, and we have 



P = {vi(x))x = 1 - {&o(x)) x - 
Therefore the rate of reproduction is 
(cri(x)cro(a; + 



(1) 



(2) 



it can be written in terms of the density-vacancy corre- 
lation function (the correlation of the of occupied and 
vacant sites) 

n m (cri(x)ao(x + 5)} x - p(l - p) 

Cw{6) = p~o-p) (3) 

Hence the rate of reproduction is 



where 



rp t {\ - Pt) 



r = l + Ci (<S=l). 



(4) 



(5) 



Obviously C\q(S = 1) is negative (a particle reduces 
the chance of its nearest neighbors to be vacant), and 
thus r < 1. In the mean- field approximation the cor- 
relation is neglected and we put r = 1. Therefore the 
mean-field equation will be 

^£l=p t (l-p t )-p t _ T (l-p t _ T ) ( t >r) (6) 

where the second term is the rate of annihilation at time 
t, equal to the rate of creation at time t—r. This equation 
is true for t > r. I assume that all existing particles at 
t — 0, gradually die during the time interval (0, r). So 
for t < t we have: 

^=p t (l-p t )-p /r (t<r). (7) 

This equation can be re-written in the integral form. 
First by integrating Eq. (7) 



pt= I p t <(l-pt<)dt' -pot/r + po (t<r) (8) 
/o 



specially for t 



Pt=r = / Pt'(l - Pt')dt'. 
Jo 



(9) 



By integrating eq. (6) and making use of cq. (9) we 
find 



Pt 



f p t ,{l-p t ,)dt' (t 

Jt-T 



(10) 



It is a definite integral with a time dependent lower- and 
upper-limit. So although the integrand is non-negative, 
p(t) may have a non-monotonic behavior. 

Finding stationary density is not possible in the dif- 
ferential equation. Setting dp t /dt — leads to nothing 
more than p t = 1 — p t - T or p t = pt- T - The former is 
irrelevant in the steady state and the latter is correct for 
every value of p. However by setting p t — pt> = p in the 
integral equation (Eq. (10) ), it turns out that 



p = 



(11) 



or 



P=l-l/r. (12) 
Thus there is a phase transition at r = r c = 1. 
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IV. OSCILLATIONS 

It is being observed that the density of particles p un- 
dergoes damped oscillations while approaching the sta- 
tionary state. The solid line in the figure 2 shows one 
such oscillatory evolution of p. Simulations are done in a 
lattice of 10000 sites with periodic boundary condition. 
The plotted curves are averaged over 100 realizations for 
t = 7. Period of oscillations is slightly greater than r. 
This kind of oscillations are present for all values of r, 
but they are weaker for smaller r. 
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FIG. 2. Density of particles versus time for 

Non-Markovian CP in real model (solid line) and in 
Mean-Field approximation (dashed line) for r = 7 and 
po = 0.75. 

The oscillatory behavior can be understood by paying 
attention to the history-dependence feature of the model. 
Since every particle dies at age r, the time evolution at 
time t is coupled to the state of the system at time t — r. 
A high creation rate at time t is equivalent to a high 
annihilation rate at time t + r. So an increase in the 
density at time t can lead to a decrease in it at some 
time later and naturally the period of the oscillations is 
of order r. 

These oscillations are also supported by the mean-field 
approximation. The mean-field equation is a delayed 
nonlinear first-order differential equation that is able to 
have oscillatory solutions which do not exist in ordinary 
first-order differential equations. Figure 2 (dashed line) 
shows one of these oscillatory solutions for r = 7. Pe- 
riod of oscillations in real model and Mean-Field are the 
same. 

The delay time in the mean-field differential equation, 
can be eliminated by a Taylor expansion of pt- r - It ba- 
sically contains derivatives up to infinite order. In fact 
here we have an infinite-order differential equation which 
is naturally able to demonstrate many complex behav- 
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iors. As in simulations, Oscillatory behavior is sensitive 
to changes in the value of r. It disappears for small 
enough values of r and it is less damped for larger r's. 



V. CRITICAL BEHAVIOR 

In this section I present the results concerning the crit- 
ical behavior of the system obtained from simulating the 
model. Simulation is made using the time-dependent 
Monte Carlo method p5| . In this method, the simula- 
tion is started with the system in a state very close to 
the absorbing state that is all the sites are vacant except 
one in origin which is occupied. The age of this single 
particle is initially set to 0. The sites are updated par- 
allel and after every time increment all existing particles 
become older by that amount. They die after growing up 
to age t. 

I measure the average population of particles N(t), av- 
eraged over all realizations, P(t), the probability of not 
entering the absorbing state up to time t, and R 2 (t) the 
mean square spreading distance. As a result of the scal- 
ing hypothesis ]l5|] , at criticality, these quantities should 
scale algebraically as 



N(t) ~ f 



P(t) ~ f 



R 2 {t)~t z 



(13) 



(14) 



(15) 



So at criticality the log-log plot of these functions should 
asymptotically become a straight line with the slope 
equal to the dynamical critical exponents. The local 
slopes for the survival probability P(t) are defined by 



5(t) 



HP(t)/p(t/b)} 

ln(b) 



(16) 



and similarly for the other exponents. I usually use 
b = 15. Away from criticality there are either upward 
or downward curvatures in the log-log plot of functions 
versus t and also in the plot of the critical exponents 
versus 1/i, depending upon the super- or sub-criticality 
of the system. By detecting the straight line from the 
curved ones, the value of r c can be evaluated with a good 
precision. Having r c , the critical exponents can also be 
found. 

Simulations are typically done up to time 400 (al- 
though many runs enter the absorbing state earlier) with 
a time increment of 0.004 for continuous-time simulation. 
Obviously this is equal to the maximum precision pos- 
sible in determining r c . Statistical quantities have been 
generally averaged over 10000 different realizations of the 
model. 
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VI. CONCLUSION 

In summary, a Non-Markovian version of the Contact 
Process is introduced. Particles have an age which de- 
termines when they are born and when they will be an- 
nihilated. Interesting oscillatory behaviors are observed 
in density of particles and the same oscillations are also 
present in the Mean-Filed approximation. Because of 
the non-Markovian property of the model there is a time 
delay in the mean-filed first-order differential equation 
which allows it to demonstrate oscillatory behaviors. 

Applying time-dependent Monte Carlo technique, crit- 
ical properties of the absorbing phase transition has been 
investigated and shown to belong to DP universality 
class. DP class is so extended to contain a non-Markovian 
model. 
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FIG. 3. Dynamical critical exponents as functions 
of 1/t. Different curves in each panel correspond to 
r = 3.06, 3.07, 3.08 respectively from top to bottom. 

Figure (3) shows the results of the dynamical simula- 
tions. In different panels local slopes as defined in eq.(16) 
for different dynamical critical exponents have been de- 
picted against 1/t. The best estimation for the critical 
age based on these graphs is t c ~ 3.07(1). For the crit- 
ical exponents I find r\ = 0.304(1), S = 0.1653(1) and 
z = 1.272(1). These critical exponents are in good agree- 
ment with those of Directed Percolation and thus the 
system belongs to the DP class. 
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